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Abstract. We generalize Galam's model of opinion spreading by introducing three competing choices. 
At each update, the population is randomly divided in groups of three agents, whose members adopt 
the opinion of the local majority. In the case of a tie, the local group adopts opinion A, B or C with 
probabilities a, (3 and (1 — a — /?) respectively. We derive the associated phase diagrams and dynamics by 
both analytical means and simulations. Polarization is always reached within very short time scales. We 
point out situations in which an initially very small minority opinion can invade the whole system. 

PACS. 82.23.Ge Dynamics of social systems - 05.65.+b Self-organized systems 



1 Introduction 

In recent years, the field of Sociophysics Jjjpl has drawn 
a growing interest from physicists, mostly theoreticians, 
with the study of opinion dynamics as one of the most ac- 
tive subjects |J,|,|l@§|Jl0lJl3. A simple model of opin- 
ion spreading with two choices was introduced some time 
ago by one of the authors In this model the system 
will eventually reach a homogeneous state called polar- 
ization by repeated debates. In this model, at each time 
step, the population is divided at random into groups of 
size m, and each group adopts the opinion of the local 
majority. This dynamics is called randomly localized with 
a local majority rule. With two choices, if m is odd, there 
are no ties and even a slight deviation from a 50/50 ini- 
tial distribution eventually leads to the polarization of the 
majority opinion ]l2] , p^[ . 

Here we investigate the richer dynamics provided by 
an m = 3 system with three possible choices, which we 
shall label by A, B, and C. This leads to the possibility 
of ties, when each member of the group has a different 
opinion |l^,^5|. We introduce the probabilities a, j3 and 
(1 — a — p) of resolving the tie in favor of A, B and C re- 
spectively. An alternative proposal for the discussion dy- 
namics has been considered in ref. ]lq| , where quantities 
analogous to a and (3 allow for a contextual interpreta- 
tion. The quantum nature of these models could reflect 
more precisely the "reduction process" by which a discus- 
sant takes up an opinion. However, in this work the focus 



a E-mail: sgekle@gmx.net 
b Associato INFN, Sezione 
pelitiOna. infn. it 
c E-mail: galam@ccr.jussieu.fr 



di Napoli. E-mail: 



lies in the reciprocal interference of opinions on three dif- 
ferent questions, while here we arc interested in a situation 
in which opinions concerning a single question can take up 
three values. 

In our model we find that polarization is always reached 
in a short time, and that it is possible to summarize the be- 
havior of the system (for given values of a and (3) by asso- 
ciating to each initial opinion distribution the value of the 
opinion which eventually prevails. We can thus draw phase 
diagrams, exhibit phase separation lines and fixed points. 
In particular it is found that in some circumstances an 
initially very small minority opinion can invade the whole 
system. While the phase boundaries are sharp in the limit 
of large populations, simulations show that they become 
blurred when the population is small, since starting from 
the same initial condition may lead to the polarization 
of different opinions. The polarization time also exhibits 
larger fluctuations, always remaining very small. 

In Section ||, after a brief review of the standard Galam 
model with two choices Q, we define the model with 
three choices and we derive the basic evolution equations. 
The fixed points and flows of these equations are dis- 
cussed in Section |[ These analytical predictions are com- 
pared against numerical simulations in Section^. Section |s| 
contains some concluding remarks and proposes some re- 
search possibilities. 



2 The model 

We first describe the randomly localized dynamics with 
local majority rule for the standard Galam model [fj^[i3| 
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with two choices, and we then introduce our present model 
with three choices. 

We consider a population of N agents, each of which 
can have opinion A or opinion B. At each time step, the 
whole population is divided into groups of size m, so that 
each agent belongs to one and only one group. Discussions 
take place, and at the end of the discussion all the mem- 
bers of the group adopt the opinion of the local majority. 
The process is iterated until a stable situation is reached. 

With m = 3, ties are ruled out and we can write out- 
right the evolution equation for the fraction p A of agents 
with opinion A ]13[: 



p A {t + 1) = P 3 A {t) + 3p A (i)(l - PA (t)). 



(1) 



This evolution exhibits the two stable fixed points at p A — 

0. pa = 1, and an unstable fixed point at p A = ^, which 
defines the phase boundary. This means that an initial 
condition slightly different form 50/50 will always lead to 
a complete polarization of the system, since the flow is 
drawn towards one of the stable fixed points. 

When there are three choices, A, B and C, there is 
the possibility of a tie in which each agent has a different 
opinion. We assume that the tie is always resolved in favor 
of one of the opinions, and we introduce the probabilities 
a, [3 and 7 = 1 — a — [3 that the winning opinion is A, 
B, or C respectively. Thus, e.g., the three members of a 
group will come out with opinion A if one of these cases 
applies: 

1. They all already have opinion A, which happens with 
probability p\ ; 

2. Two of them have opinion A, which happens with 
probability 3p A (I — p A )\ 

3. There is a tie resolved in favor of A, which happens 
with probability 6ap A pB(l — p A — Pb)- 

A similar analysis can be made for pb- We consider the 
large population limit, in which sampling fluctuations can 
be neglected. We thus obtain the evolution equations for 
the fractions p A , pb- 

p A (t+l)=p 3 A {t) + 3p 2 A (t)(l-p A (t)) 

+ 6a PA (t)p B (t)(l ~ PA (t) - PB (t)); (2) 
p B (t+l)=p B (t)+3pi(t)(l-p B (t)) 

+ 6(3p A (t)pB(t)(l-p A (t)- P B(t)). (3) 

The state of the system can be conveniently repre- 
sented by a point in an equilateral triangle. Let A, B 
and C be the vertices of such a triangle, then the state 
defined by the probabilities p A , pb and pc — f — p A — Pb 
is represented by the point 



p = p A A + p h B + p c C 
= C+paCA + pbCB, 



(4) 



where CA and CB are the vectors leading from C to A 
and B respectively. The vertices of the triangle represent 
the polarized states. Inspection of eqs. (pP) (and com- 
mon sense) shows that these states are fixed and stable. 
Each vertex thus commands a nonempty attraction basin. 



Indeed, one can see that if, e.g., \ < p A (t) < 1, one has 
p A (t + l) > p A (t). Thus the attraction basin of A contains 
at least the region p A > i . 



3 Fixed points and flow diagrams 

Equations (0J|) exhibit the following fixed points: 



I. 



The trivial fixed points A, B and C, corresponding to 
the polarized states; 

Since, when one of the three fractions p A , pb or pc 
vanishes, the evolution equations reduce to the one de- 
scribing the two-choice Galam model (eq. (Q), there are 
the corresponding fixed points D (p A — 5, pb = 5), E 
(pb = 5, PC = 5), and F (p c = |, p A = ±) which lie 
on the triangle sides. 

Moreover a seventh fixed point G, in which all three 
fractions are nonzero, may appear in the DEF triangle. 
For this point to appear, it is necessary that none of 
the probabilities a, (3 and 7 = I — a — f3 is larger than 
|. Thus the region in the (a, j3) plane in which this G 
point appears is the one dashed in fig. [l| 



2/3 



1/3 



Fig. 1. The region in the (a, (3) plane in which the inner fixed 
point G appears is dashed. 



The stability of the fixed points can as usual be dis- 
cussed by linearizing the evolution equations around them. 
If we denote by J the Jacobian matrix of the transforma- 
tions (|^. |) , a fixed point is stable if all the eigenvalues of J 
are smaller than I in absolute value. One can easily check 
that: 

1. The three trivial fixed points A, B and C have all eigen- 
values equal to 0; 

2. The three side points D, E and F, have one eigen- 
value equal to |, corresponding to the eigenvector ly- 
ing along the side of the triangle, whereas the other 
one is given by 

^(1 - a - (3) for D, for E, and |^ for F. (5) 
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We see that the direction towards the interior of the 
triangle becomes unstable if the probability parameter 
associated with the opposite vertex is larger than §. 
In this situation, the inner fixed point does not exist. 
3. The inner fixed point, if it exists, has both eigenvalues 
larger than 1 in modulus. 

The basic structure of the flow is thus represented in fig. ||. 
The first applies to the case in which there is the inner 
fixed point, G, (Case I) and the other to the case in which 
it does not exist (Case II). It is clear from this figure that, 
in Case I, the attraction basins of the stable fixed points 
meet two by two at the fixed points D, E and F, and 
all three in the inner point G. In Case II, on the other 
hand, the attraction basin of the point corresponding to 
the largest favorable probability (a for A, etc.) will sepa- 
rate the basins of the other fixed points at the side fixed 
point situated midway between the two. 

The phase boundaries can be conveniently located by 
iterating the inverse of eqs. (||J^), starting close to the side 
fixed points D, E, or F. The iterations converge to the in- 
ner fixed point G, if it exists, or to the other relevant side 
fixed point. As an example, we show in fig. ^ the phase 
diagram obtained for a = = ^ and | respectively. It 
is clear that the attraction basins of the side fixed points 
reduce to stretches of the phase boundaries: the one con- 
necting each point to G, in Case I, or to the unstable side 
fixedpoint (like D in fig. |L right), in Case II. In figures [i] 
and |j we show respectively the flow and the phase dia- 
grams for a — 1/2, f3 = 1/3 (Case I) and for a = 2/9, 
(3 — 1/9. We can conclude that in Case II, close to the 
unstable side fixed point, there is the possibility that the 
"preferred" opinion (the one with the largest probabil- 
ity of resolving a tie) can invade the system even if it is 
initially supported by a very small minority. The width 
of the sector in the p = (j>a,Pb) plane where this can 
happen rapidly vanishes as p approaches the side of the 
triangle, since the two boundary lines are tangent to each 
other (unless a — [3 — 0). The sector, however, opens up 
more and more rapidly as the third opinion becomes more 
and more preferred. 

If one of the three probabilities (e.g., (3) vanishes, the 
straight line pb = \, connecting the side fixed points D 
and E on either side of the corresponding vertex B, is 
conserved by the flow and defines the boundary of the at- 
traction basin for it. This means that the "unpreferred" 
opinion can invade only if it is initially shared by a major- 
ity. For a < | the fixed point E is stable along this line. 
For i < a < | the inner fixed point G appears on the line, 
and both side fixed points D and E are partially stable. 
Finally, if a > |, the fixed point D becomes partially sta- 
ble, while E becomes totally unstable. The different cases 
are represented in fig. ||. 

4 Simulations 

4.1 Phase diagram 

We now consider the behavior of systems with finite values 
of N. Equations (§,||) provide a good description of its 



behavior, provided one decrees that complete polarization 
has taken place whenever 

Pi<^, ie{A,B,C}. (6) 

In this case, indeed, every single agent supports the win- 
ning opinion. For TV large enough, the attraction basins we 
have identified on the basis of the iteration equations hold 
with practical certainty, except very close to the bound- 
aries. When N gets smaller, one can see some uncertainty 
building up along the boundaries, as shown in fig. |?]. Here 
we assign a color to each of the vertices of the triangle 
(black for A, light grey for B, and dark grey for C) and we 
assign to each of the initial conditions (a point (pa>Pb) 
in the triangle) the color of the vertex reached most often 
in 50 trials. The diagram is evaluated on approximately 
5000 different initial conditions with different values of N: 
N = 10 4 , N = 2500, and N = 100. It is clear that while 
the overall shape of the phase diagram remains unchanged, 
the phase boundaries become more and more blurred as 
N decreases. 

To get a more precise measure of the difference be- 
tween the results of simulations and the predictions of the 
deterministic equations, we compared the phase diagrams 
one by one for all 5100 initial distributions and counted 
the percentage of mismatches. Results are shown in ta- 
ble 1. One can see how the quality of agreement diminishes 
with smaller system sizes, but can be improved again by 
increasing the number of runs over which the statistical 
average is taken. This improvement only works up to a 
certain level, though, the error margins for 5000 runs in a 
system with N = 100 stay almost exactly the same as for 
1000 runs and are therefore not shown here. As one can 
see in fig. fj], the mismatches occur only around the phase 
separation lines. As soon as the initial distribution is set 
a small distance apart from these lines, the percentage of 
mismatch goes towards 0. As the error even for small sys- 
tems with N = 100 and only 50 runs is maximal 2.4%, we 
conclude that the deterministic eqs. (|],[|) work very well 
for phase diagrams down to this size. 



4.2 Polarization time 

It is also interesting to discuss the polarization time T p , 
i.e., the number of discussion cycles needed to reach com- 
plete polarization. It depends on the initial condition as 
well as a and (3. We list in table polarization times for 
a few interesting configurations. In fig. || the polarization 
times are displayed as a function of the initial condition 
over the whole triangle. One can see that the polarization 
time becomes longer as the initial distributions approaches 
the phase boundaries. The scaling of T p with system size 
N, using the deterministic evolution equations, is shown 
in fig. H for two configurations, one with an initial distri- 
bution close to the triple point G with a = [3 = starting 
from pa = 0.51 and p-g, = 0.49, the other far from the 
triple point with a = [3 = k starting from pa = Pb = 0.1. 
One sees that, even if N varies from 25 to 10 s , T p remains 
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Fig. 2. Fixed points and basic structure of the flow. Left: Case I, in which (a, /3) lies inside of the dashed region in fig. |l|. Right: 
Case II, in which (a, j3) lies outside that region. In the case shown, 7 = ! — a — /3> |. 





Fig. 3. Phase diagrams for the dynamics: Left: a = /3 = 4 (Case I). Right: a = f3 = | (Case II). 



N = 10 4 
(50 runs) 



N = 2500 
(50 runs) 



N = 100 
(50 runs) 



N = 100 
(100 runs) 



N = 100 
(1000 runs) 



a=± jy 3=§ 
a = 0, (i = 
a= |, /3 = 
a=| )( S = 
a= |, /3 = 



0.7% 
0.1% 
0.5% 
0.1% 
0.1% 



1.0% 
2.0% 
1.5% 
1.1% 
1.2% 



1.8% 
2.4% 
2.2% 
2.0% 
1.8% 



1.5% 
2.1% 
1.9% 
1.7% 
1.6% 



1.0% 
2.1% 
1.6% 
1.3% 
1.4% 



Table 1. Percentage of mismatching points between probability calculations and simulations for different system sizes. 



virtually constant, varying only between 11 and 14 in the 
first and between 2 and 5 in the second case. 

The scaling of the simulated polarization time for the 
same configurations as above is shown in fig. [h]. The pic- 
ture now looks quite different than for the phase diagrams. 
Regarding T p only for systems with about 1000 agents the 
agreement is still acceptable, although the error bars are 
already quite large. Overall the agreement for the system 
starting far from the triple point is much better than for 
the system starting close to it. 



Looking at the relative error shown in fig. Ill], one must 
suggest that for small systems with TV < 1000 a predic- 
tion of the polarization time via the evolution equations 
(^],[5]) is unreliable. However, predictions about the winning 
opinion remain fairly reliable as described in the previous 
section. 
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Fig. 4. Flow, fixed points and attraction basins for a = 4/9, /3 = 1/9 (Case I). Left: Flow and fixed points. Right: Boundaries 
of the attraction basin. 





Fig. 5. Flow, fixed points and attraction basins for a — 2/9, /? = 1/9 (Case II). Left: Flow and fixed points. Right: Boundaries 
of the attraction basin. 






Fig. 6. Phase diagrams with (5 
Right: a = §. 



0. The straight line DE (p B 



is conserved by the flow. Left: a 



j_ 

10 ■ 



Center: a 



_6_. 

12 ' 



5 Concluding remarks 

In this paper we presented an extension of Galam's model 
for a two-choice system Jl3| to a system in which each 
agent can choose between three possible opinions. All re- 
search has been done for discussion groups of size 3 and 
in a randomly-localized model. To resolve ties, two new 
parameters a and (3 measuring "argument strength" were 
introduced. Opinion distributions could be conveniently 
visualized in an equilateral triangle. 



We showed that, independent of the parameter values, 
the system will almost always reach a complete polariza- 
tion after a small number of discussion cycles. These po- 
larization times are virtually independent of system size. 
We have drawn the phase diagrams, associating to each 
initial condition (pai.Pb) the opinion (A, B, or C) which 
eventually wins. Finally, a number of simulations has been 
conducted to show that the deterministic evolution equa- 
tions used to obtain the above results are reliable down 
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Fig. 7. Simulated phase diagrams for (from left to right) a — f3 = j; a = (3 = 0; a = \ and (3 = 0; a = | and /3 = 0; a = | 
and /9 = 0. From top to bottom system sizes of N = 10 4 , N = 2500 and N = 100. The color of a point corresponds to the 
vertex reached most often in 50 runs, starting from that initial condition. One can see how the boundary lines become blurred 
for small systems. 













time T p 


time T p 


a 


PA 


pB 


winner 


(probabilities) 


(simulated) 








1 





A 





± 








0.51 


0.49 


A 


13 


13.1 ± 1.4 








0.49 


0.49 


C 


13 


13.8 ± 1.1 








0.30 


0.30 


C 


5 


5.4 ± 0.5 








0.1 


0.1 


c 


3 


3.1 ± 0.3 


1/3 


1/3 


0.1 


0.1 


c 


4 


4.0 ± 0.3 


1/2 





0.1 


0.1 


c 


4 


4.0 ± 0.2 


1/2 





0.3 


0.51 


B 


13 


13.2 ± 1.4 



Table 2. Polarization time T p for a few interesting configurations. Results obtained by the evolution equations (|2j|3j) and by 
simulations for a system of = 10 4 . Simulated results averaged over 100 runs. 



to systems with approximately 100 agents, as far as the 
phase diagram is concerned, and to N ~ 1000 for the po- 
larization times. 

Similar phase diagrams are obtained with a serial up- 
date dynamics, in which, at each microscopic time step, 
three different agents are picked up at random and their 
opinions updated according to the local majority rule de- 
scribed above. One counts a time step when each individ- 
ual has taken part to three discussions on average. The 
polarization times with this dynamics are longer. Indeed, 
while in the parallel dynamics discussed so far the ap- 
proach to the stable fixed points is ultrafast (the deviation 
gets squared at each time step) it is only exponential in 
the serial update dynamics. Thus 



Vh7]v, 

IniV, 



with parallel dynamics; 
with serial dynamics. 



(7) 



Our results show that dealing with three opinions is 
much more complex than with two. From the complexity 
of the various flows, it is seen that initial conditions are 
not enough to predict any outcomes since the values of the 
parameters a and (3 are instrumental to determine the 
highly nonlinear road to success or failure in the rather 
quick competition in the opinion forming process. What 
appears as a successful strategy for one opinion in the first 
steps may turn out to be a total disaster at the end. At 
this stage of our investigation the following clues can be 
outlined. 

1. It seems that we have here some hint on why the var- 
ious attempts to create a third "way" in democratic 
countries have always failed. Starting from an A-B two- 
party situation two possibilities exist for the creation 
of a third one C. Either it is created on its own against 
A and B or it is "helped" by one of the two big forces, 
say A. In the first case we have a = (3 — ^ making C 
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Y \\ V 
Y V 

Fie. 8. Polarization time T p in a system with N — 10 4 for different parameters obtained by the deterministic evolution equations 
(^|j). Darker colors mean longer time. Note the peaks around the phase boundaries. Above: left: a = j3 — |; center: a — j3 — 0; 
right: a = ~, /3 = 0. Below: left: a = ~, /3 = 0; right: a = §, /? = 0. 



a = /3 = 
a = (3 = 1/3 



a = /3 = 
a = /? = 1/3 A 



10000 100000 1e+06 

System size 



1000 
System size 



Fig. 9. Polarization times T p obtained by iteration of eqs. 
for different system sizes. The system starts fom pa = 0.51, 
Pb = 0.49 for the case a = f3 = (squares) and from pa = 0.1, 
Pb =0.1 for a = j3 = 1/3 (triangles). The same convention ap- 
plies to figs. ^| and [H]. The system starting far from the triple 
point G (lower curve) converges a lot faster on a polarization. 
Even over a wide range of system sizes T p remains virtually 
constant. 



much further from the boundaries since by its nature 
it is a new comer. Therefore it will rather quickly dis- 
appear. In the second case we have a = (3 = since at 
a tie, A goes for C. In this case C will start growing 
slowly at the expense of both B and A. To stop it from 
spreading all over the unique solution is to have A and 
B to make a coalition against C with a — (3 = |. 
Focusing on the polarization effect with one opinion 
spreading in two opinion competitions, it seems that 
the best strategy is not to reinforce its own side but 



Fig. 10. Polarization times T p obtained through simulations 
for different system sizes. As expected, the error bars become 
smaller with increasing system size. Overall, error bars for the 
system starting far from the triple point G are much smaller 
than for the other system. Results averaged over 100 runs. 



instead to create a third opinion which cooperates with 
either one of the others. It makes the new very minority 
opinion to eventually win the competition. 

The above ideas are of course qualitative and need to be 
deepened. This paper is just a beginning outlining some 
basic features of three-opinion systems. There is much 
room for further research, as for example the introduction 
of a geographical distribution and the influence of neigh- 
borhood relations Q. Another interesting topic would be 
the behavior of discussion groups with more than three 
members 14 . To introduce some dependence of the pa- 
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10 100 1000 10000 100000 

System size 

Fig. 11. Relative error of the simulated polarization time T p . 
For systems with N < 1000 the error becomes very large and 
predictions seem unreliable. 



16. D. Aerts and S. Aerts, Foundations of Science 1, (1995) 85, 
reprinted in Topics in the Foundation of Statistics, edited by 
B. Van Fraassen (Kluwer, Dordrecht, 1997). 



rameters a and j3 on the values of the respective support 
for A and B would be also of a great interest. 
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